Observable primordial vector modes 
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Primordial vector modes describe vortical fluid perturbations in the early universe. A regular 
solution exists with constant non-zero radiation vorticities on super-horizon scales. Baryons are 
tightly coupled to the photons, and the baryon velocity only decays by an order unity factor by 
recombination, leading to an observable CMB anisotropy signature via the Doppler effect. There is 
also a large B-mode CMB polarization signal, with significant power on scales larger than I ~ 2000. 
This B-mode signature is distinct from that expected from tensor modes or gravitational lensing, 
and makes a primordial vector to scalar mode power ratio ~ 10~® detectable. Future observations 
aimed at detecting large scale B-modes from gravitational waves will also be sensitive to regular 
vector modes at around this level. 
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Observations of the eosmic microwave background 
(CMB) show that the primordial perturbation was al- 
most certainly dominated by adiabatic scalar (density) 
modes. However it is well known that there are sev- 
eral possible scalar isocurvature modes T| that could be 
present at some level. In the presence of a primordial 
magnetic field, there is also an observable vector mode 
perturbation sourced by the anisotropic stress of 
the magnetic field. Other sources such as topological de- 
fects can also source vector modes. Here we concentrate 
on the rarely-considered regular primordial (unsourced) 
vector modes, which are non-decaying solutions of the 
perturbation equations in the presence of free streaming 
neutrinos |^] . We show that a very small primordial reg- 
ular vector mode amplitude could be observable. 

In the absence of an initial large scale radiation vortic- 
ity the vector modes remain in a decaying mode and have 
essentially no observational signature. They are therefore 
not predicted to be present at any significant level in in- 
flation or other simple models. However there is a regu- 
lar mode with a non-zero initial photon vorticity, having 
equal and opposite initial photon and neutrino angular 
momenta such that the total large scale angular momen- 
tum is zero. This is the vector analogue of the scalar neu- 
trino isocurvature velocity mode discussed in Ref. Q , and 
constitutes a valid possible component of the general pri- 
mordial perturbation. These velocity modes would have 
to be excited after neutrino decoupling and are hence 
difficult to produce and somewhat contrived. But they 
remain a logical possibility that can be constrained by 
observation, and if observed would be a powerful way 
to rule out most theoretical models (for constraints on 
the scalar mode see e.g. Ref. and references therein). 
The vector mode can be detected at very small ampli- 
tudes and distinguished from the various scalar modes 
because of its distinct non-zero B-mode (curl-like) CMB 
polarization signal that is absent with only linear scalar 
modes. 

As we show, the vector _B-mode signature is quite dif- 
ferent from that expected from weak lensing or primor- 



dial tensor modes. On large scales the spectrum is similar 
to that from tensors, so observations aimed at detecting 
the i3-modes from primordial tensors will also be sensi- 
tive to the large scale part of the vector power spectrum, 
but they can easily be distinguished by the vector mode 
power on smaller scales. The physical difference between 
the spectra is that tensor modes rapidly decay as soon as 
they come inside the horizon, whereas the vortical modes 
are nearly constant during radiation domination, decay- 
ing only on small scales though damping towards the end 
of tight coupling. 

Even more contrived regular modes exist with non-zero 
primordial neutrino octopole (or higher) 0, j however 
these have a much weaker observational signature and 
are not considered further here. 



A. Covariant Equations 

We consider linear perturbations in a flat FRW uni- 
verse evolving according to general relativity with a cos- 
mological constant, neglect any velocity dispersion of 
the dark matter and baryon components, and approxi- 
mate the neutrinos as massless. Perturbations can be 
described covariantly in terms of a 3-1-1 decomposition 
with respect to some choice of observer velocity Ua (we 
use natural units, and the signature where Uau"- = 1), 
following Refs. 01 13, 0- Projected spatial derivatives 
orthogonal to Ua can be used to quantify perturbations 
to scalar quantities, for example the pressure perturba- 
tion can be described in terms of D^p where the spatial 
derivative is 



Da = Va - UaU Vfc. 



(1) 



Conservation of total stress-energy V°Tab = implies 
an evolution equation for the total heat flux qa 



Qa + ^Oqa + {P + p)Aa - BaP + D^TTafc = 0. 



(2) 
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where p is the energy density, qa = u^^i,qa, O = Vita 
is three times the Hubble expansion, Aa = UbV^Ua is 
the acceleration, and Tr^b = T^afc) is the total anisotropic 
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stress. Angle brackets around indices denote the pro- 
jected symmetric trace-tree part (orthogonal to Ua). 

We define the vorticity vector il^ = curl Ua where for 
a general tensor 



curl Xr, 



Vbcdiay^D'X^,^,,,,^^ (3) 



and round brackets denote symmetrization. It has the 
evolution equation 



(4) 



and is transverse DTia = 0. Remaining quantities we 
shall need are the 'electric' Eab and 'magnetic' Hat parts 
of the Weyl tensor Cabcd 



Eab = CacbdU'^u'^ 



1 , , 

Hab = -^VacdfCbe U'^U^ (5) 



(which are frame invariant) and the shear aab = D^^jWf,^. 
The Einstein equation and the Bianchi identity give the 
constraint equations 

1 2 

D^CTah - - cm-m,, - gD^e - Kqb = 

^"Eab - n (^qb + ^Dbp + \^'"^ab \ = 



1 



+ P)^b + curl qb] = 



Ha 



curl Un 



(a^^b) 



(6) 



and the evolution equations 



Eab 
Hab 



-QEab 
QHab 



1 

= —Eab — -^K-Kab 
= curl Hab + ^ 



= - curl Eab - 2 curl T^ab ■ 



(7) 



D^aSf,) for some first 



We use natural units where c = 1 and define n = SttG. 

A vector like Aa may be split into a scalar part A'^a'^ 
and a vector part A'}^ where Aa = ^i"'' -I- A^'p , aI"^ = 
Dav4 for some first order scalar A and the vector part is 
solenoidal D'^A^a^ = 0. This extends to a tensor where 
the vector part is given by cr^jj'' 
order solenoidal vector S;,. 

We can choose Ua to simplify the analysis. At lin- 
ear order one can always write Ua = -I- Wq, where 
is hypersurface orthogonal and Va is first order, so 
curlita = curlua. For a zero order scalar quantity X 
it follows that DqAT — D^AT — VaX. For vector modes 
{Q^X)'^^^ — 0, and it is convenient to choose the frame 
Ua to be hypersurface orthogonal so that curl = and 
hence (DaA)^^) — 0, where the bar denotes evaluation in 
the zero vorticity frame. From its propagation equation, 
vanishing of the vorticity also implies that A'^p = 0, so 
the zero vorticity frame coincides with the synchronous 



gauge. The CDM velocity is also zero in this frame mod- 
ulo a mode which decays as 1/S' where S is the scale 
factor. 

It is convenient to expand the vector components 
in terms of transverse eigenfunctions of the zero order 
Laplacian, where S^D^Q^ — k^Qa and ± denotes 
the parity. A rank-£ tensor may be expanded in terms of 
rank-^ eigenfunctions defined by 



Q 



D 



(ai 



which satisfy 



k {P - 1) 



n^'O* — — 

VA,_ia, S 1(21-1)^'''-' 

1 k 

curlQ± = Js^X- 
Harmonic coefficients are defined by 



(8) 

(9) 
(10) 



k 



(11) 



where the k and ± dependence of the harmonic coeffi- 
cients is suppressed and qi = {pi -\- Pi)vi for each fluid 
component, where Vi is the velocity, and the total heat 
flux is given by q = Ylili- The sum is over k and the ± 
parities. We write the baryon velocity simply as v. 

The equations for the harmonic coefficients in the zero 
vorticity frame reduce to 



k{a + 2na) 



2KS^q = k^a 



(12) 



where the dash denotes a derivative with respect to con- 
formal time T], and H — 5*9/3 is the comoving Hub- 
ble parameter. The combination v + a (the Newtonian 
gauge velocity) is frame invariant, as are a — a + and 
V — V ~ n. By choosing to consider the zero vortic- 
ity frame we have simply expedited the derivation of the 
above frame invariant equations. Other papers use the 
Newtonian gauge 10], in which a is the vorticity. 

The evolution equation for the shear has the solution 



-1 f , KS^n 



(13) 



In the absence of anisotropic stress it therefore decays as 
1/5^. However after neutrino decoupling the neutrinos 
will supply an anisotropic source, and solution of this 
equation requires a consistent solution for the neutrino 
evolution. 

The baryon velocity is coupled to the photon velocity 
via Thomson scattering 



Tiv — ——SUeCFT ( -V — Ii 

Pb 



(14) 
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where Ji = 4z;^/3, rig is the electron number density and 
ctt is the Thomson scattering cross-section. 

The photon multipole equations for vectors become 



I',+k 



I 



21 + 1 



(^ + 2) 



4 _ 2 , ^ _ 



where C, = 3/2/4 — 9E2/2 is a source from the photon 
anisotropic stress and i?-mode polarization, and — 
n^. In general /; is an angular moment of the frac- 
tional photon density distribution, four times the frac- 
tional temperature anisotropy. The neutrino multipole 
equations are analogous but without the Thomson scat- 
tering terms. The solution is 



IiiVo) 



no 



drje 



Suearv^iix) 
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fccr 



dx 



(16) 



where ^;(a;) = lji{x)/x, ji{x) is a spherical Bessel func- 
tion, X = ^(^70 ~ v) ^-iid T is the optical depth. In the 
approximation that the visibility SneCTTS.~^ is a delta 
function at last scattering rj = rf this becomes 



cm 

4 dx 



+ 2 driH'-^i. (17) 

r]* n* 



The anisotropy therefore comes predominantly from the 
Newtonian gauge baryon velocity at last scattering, plus 
an integrated Sachs- Wolfe (ISW) term from the evolution 
of the magnetic Weyl tensor H — a/2 along the line of 
sight. 

The vector polarization multipole equations prj be- 
come 



Br 



(Z + 3)(Z + 2)Z(/-1) , 

'+1 21 + 



(? + 1)3(2/ + 1) 
2 



-kE, 



+ 



{l + 3)il + 2)l{l~l) 
{I + If (21 + 1) 

I(7TT)^^'"^° 
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-kBf , 
1 



(18) 



where Ei and Bi describe moments of the E (gradient- 
like) and B (curl-like) polarization. These equations have 
solutions 



I - 1 

i-i 



dr]SneiJTe 



d^iix) , 2*i(x) 



dx 



c 



I - 1 r° 

j-^ j dTjSneare-^^iixK- 



(19) 



Signs of El and Bi here follow the conventions of CMB- 
FAST [13 and CAMB Q^. 
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FIG. 1: Evolution of the vector perturbations with wavenum- 
ber k = 0.02Mpc"^ (solid lines) and k = 0.2Mpc"^ (dashed 
lines). Thick lines are the velocities of the baryons (top), 
photons (until decoupling same as the baryon velocity) and 
neutrinos (bottom). Thin lines are a (with ao = 1) and the 
photon anisotropic stress. The baryon velocity evolves inde- 
pendently of wavenumber on large scales, but is damped on 
small scales. 



B. Solutions 

At early times the baryons and photons are tightly 
coupled, the opacity = Sn^aT is large. This means 
f-y ~ t;, and we can do an expansion in Tc that is valid 
for e = max(A;Tc, TYtc) <C 1. To lowest order 



RHv 

1 + R 



+ 0(T-e) 



where R = Zp^/'ip^. The solution is 



Wo 



1-1- i? 



(20) 



(21) 



where vq is the initial value. Hence if vq 7^ 0, by de- 
coupling V has only decayed an order unity factor de- 
pending on the matter and radiation density at the time. 
On smaller scales where kxc = 0{\) before decoupling 
the perturbations are damped by photon diffusion, giv- 
ing a characteristic fall off in perturbation power on small 
scales. 

We now perform a general series expansion in confor- 
mal time for the above equations in the early radiation 
dominated era to identify the regular primordial modes. 
We define uj = QrriHo / V^R, where flR = + and 
Ti.0 and fli are the Hubble parameter and densities (in 
units of the critical density) today. The Fricdmann 
equation gives 



S = 



(22) 



4 




FIG. 2: Typical CMB temperature (top solid), polarization 
EE (bottom solid), BB (dashed thick) and cross-correlation 
TE (dotted; absolute value) power spectra for regular vector 
modes assuming a primordial vector to scalar power ratio ~ 
10~^ and scale invariant vector mode spectrum P^. The other 
dashed lines show the B-mode spectrum from weak lensing 
(peaking at £ ^ 1000), and primordial tensors with initial 
power ratio ^ 10~^ (peaking at ~ 100). 



Defining the ratios = ^u/^Ri Rj = ^^7/^^^, Rb = 
fib I , and keeping lowest order terms the regular solu- 
tion (with zero initial anisotropies for Z > 2) is 



(70 I 1 



CTq- 



R^ 
aoAR, 
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YaR^ + 15 
5/1 3Rb 
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(23) 
(24) 
(25) 
(26) 



where we have neglected small contributions from the 
scattering-suppressed photon anisotropic stress. This 
regular mode is the vector analogue of the neutrino ve- 
locity isocurvature mode discussed in Ref . . The shear 
a is initially constant on super-horizon scales, supported 
by the growing anisotropic stress of the neutrinos. On 
sub-horizon scales in radiation domination it decays as 
the neutrino anisotropic stress starts to oscillate rather 
than grow. 

The photon and neutrino vorticities are constant on 
super-horizon scales during radiation domination. This is 
consistent with angular momentum conservation because 
of the energy redshift. The photon vorticity is tightly 
coupled to the baryons, so both are initially nearly con- 
stant, with some decay due to drag from the baryons 
through matter radiation equality. On super-horizon 
scales there is only an order unity decay, so a signif- 
icant large scale photon quadrupole will be present at 
low redshift to source a significant additional large scale 



polarization signal from reionization. The evolution is 
illustrated in Fig. ^ 

On large scales the early ISW contribution is about 
20% as tr decays as the matter becomes more dominant. 
On scales sub-horizon at recombination there is no ISW 
contribution as a has already decayed. We neglect the 
effect of magnetic field generation by the photon-baryon 
vorticity Q. 



C. Observations 

We now compute the observable CMB anisotropy sig- 
nal. We define the dimensionless first order transverse 
vector Oa such that a'^^jj = D(^aO'b) , and quantify the pri- 
mordial vector modes by their power spectrum defined 
so that 



{\^a\ 



dink 



(27) 



The corresponding expressions for the CMB temperature 
and polarization power spectra are derived in Q . 

To account for the small scale damping effect accu- 
rately, as well as a detailed treatment of recombina- 
tion and reionization, we compute sample CMB power 
spectra numerically by a straightforward modification of 
CAMB^ The CMB power spectra (Q) depend on Pg. 
For a scale invariant spectrum, the temperature Ci has 
a broad peak around £ ~ 50, as shown in Fig. |2 The 
polarization power spectra peak at around £ ~ 500, with 
the _B-mode dominating in accordance with Ref. |l(]|. 

The large scale reionization signal is rather similar to 
that expected from tensor modes, and thus experiments 
aimed at detecting this tensor signal will also be sensitive 
to vector modes. Incomplete sky coverage only decreases 
the sensitivity by an order unity factor due to E-B mode 
mixing [T^ [l5l | even on the largest scales. From Fig. [21 
we see that the large scale B-modes are more sensitive 
to vector power by a factor of about 100, thus sensitive 
observations of tensor modes will also be good probes of 
regular vector modes. To distinguish the two one just 
needs to measure the spectrum at ^ > 100 where the 
tensor power falls but the vector power continues to grow. 

The dominant confusion on small scales is likely to be 
from weak lensing of the scalar modes, which peaks on 
similar scales. There are about 10^ observable modes, 
so one can ideally expect to detect a vector contribu- 
tion ^ 1/1000 of the power of the lensing signal. Since 
they are of comparable power for a scale invariant pri- 
mordial power spectrum ratio P^/P^ of ~ 10^"^ [P^ 
is the power in the comoving curvature perturbation), 
this implies that vector modes with only 10~^ of the 
scalar power may be detectable irrespective of the ten- 
sor mode amplitude. Since the lensing signal is non- 
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Gaussian, and in the absence of vector modes is partially 
subtractable 0, 0, 0| , the in-principle limit is proba- 
bly much lower, though this depends on the spectrum of 
the vector modes. The ultimate limit may be around the 
level where there should be a sourced vector mode signal 
from second order effects fl^ - 

Primordial magnetic fields source a B-mode spectrum 
similar to that from primordial vector modes 2]. How- 
ever the perturbations are expected to be highly non- 
Gaussian for magnetic fields and hence easily distinguish- 
able from primordial vector modes if they are approx- 
imately Gaussian, at least until the lensing confusion 
limit. Magnetic fields also provide a constant source 
which partly compensate the damping, so there is more 
magnetic field vector mode power on very small scales. 
The detailed signature of magnetic fields in the CMB is 
discussed in Ref 0, including the additional large scale 
signature from tensor modes. 

Topological defects can also source similar B-mode 
spectra (20j, though again the spectrum is expected to 
be non-Gaussian, and (at least for strings) there is more 
power on very small scales due to the continuous sourcing 
of the vector modes. 



bility that tiny primordial amplitudes can be constrained 
from future high-sensitivity CMB polarization i?-mode 
observations. Any signature of vector modes would be 
powerful evidence against simple inflationary models. 
The Planck^ satellite should be able to detect the B- 
mode signature from primordial vector modes at the 10"'^ 
level, and distinguish them from tensor modes by the 
presence of small scale power. A full Bayesian joint anal- 
ysis of all the CMB power spectra should be straightfor- 
ward using MCMC techniques, and may give better con- 
straints that suggested here. Separating a vector mode 
signal at the 10~^ level from that generated by lensing of 
scalar modes would be a serious challenge for the future. 
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